Coupled Microwave Billiards as a Model for Symmetry Breaking 
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Two superconducting microwave billiards have been electromagnetically coupled in a variable 
way. The spectrum of the entire system has been measured and the spectral statistics analyzed as 
a function of the coupling strength. It is shown that the results can be understood in terms of a 
random matrix model of quantum mechanical symmetry breaking - as e.g. the violation of parity 
or isospin in nuclear physics. 

PACS number(s): 05.45.+b, 11.30.Er 



Both, classical and quantum mechanical chaos can be 
studied with the help of billiards, see e.g. the recent re- 
view Q| . Quantum mechanical billiards are readily simu- 
lated by sufficiently flat microwave resonators since 
the Schrodinger and the Helmholtz equations are equiv- 
alent in two dimensions. Here, we present a study of a 
system consisting of two coupled resonators. This sim- 
ulates the breaking of a symmetry, see below. The res- 
onators were made superconducting. This allowed us to 
study the transition from the uncoupled case into the 
regime of weak coupling very precisely. The positions of 
the resonances were determined with a precision of 10~ 7 . 

Symmetry breaking in chaotic systems has been in- 



tensely investigated. Im 
violation in heavy nuclei! 
metries were studied in Ifi 



jressive is the study of parity 
3fl. Atomic and molecular sym- 
M. Another example from nu- 
clear physics is isospin mixing, see e.g. [p|jic|l. In jt(|, the 
complete spectrum of the nucleus 26 Al at low excitation 
energy was established. The analysis of this spectrum 
in terms of the so called Gaussian Orthogonal Ensemble 
(GOE) of random matrices showed that the level statis- 
tics was intermediate between a 2-GOE and a 1-GOE 
behavior |fTl]j . By this, we refer to the following model. 
Each level in the spectrum of 26 Al can be characterized 
by isospin or 1. In the absence of mixing, the spectrum 
of the states of each isospin (i.e. of each symmetry class) 
has the statistical properties of the eigenvalues of matri- 
ces belonging to the GOE. The superposition of the two 
spectra displays a 2-GOE behavior. It is described by 
the first term of the Hamiltonian 
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This is a special case of the model of Ref. m. The off- 
diagonal matrix couples the classes. The random ele- 
ments in the GOEs and in V , all have the same rms value 
v so that a = 1 makes 7i as a whole to be a GOE ma- 
trix. The resulting spectrum displays 1-GOE behavior. 
For the observables studied below, the 1-GOE behavior 



is actually reached already if av/D is w 1. Here, D is the 
mean level distance of 7i. For simplicity, we set v = 1 in 
the sequel. This makes D dimensionless and the param- 
eter governing the level statistics is then a/D. - In the 
example of 26 Al, this parameter was determined from the 
level statistics; whence the mean square Coulomb matrix 
element that breaks isospin was derived. The present 
experiment tests the model (|l|) with a large number of 
levels and very clean spectra in a situation, where the pa- 
rameter a/D controlling the symmetry breaking can be 
varied. Alternative models for coupled chaotic systems 
are e.g. given in 0Q, 

7i=l 7 z =l-8 
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FIG. 1. Shapes and locations of the antennas of the two 

coupled Bunimovich stadium billiards. 

In the present experiment each of the two symmetry 
classes consists of the eigenstates of a (quarter of a) sta- 
dium billiard, see Fig. [l]. The radius of the quarter circle 
was in both cases r = 0.2 m. The ratios 7 between the 
length of the rectangular part and r were 71 = 1 and 
72 = 1.8 for the two billiards, respectively. The measure- 
ment was restricted to frequencies below 16 GHz where 
both resonators are two-dimensional and display 608 and 
883 resonances, respectively. For the variable supercon- 
ducting coupling, the two resonators were put on top of 
each other and holes, 4 mm in diameter, were drilled 
through the 2 mm thick walls of both resonators (see 
Fig. ||). A niobium washer ensured sufficient electrical 
contact between the resonators. Coupling was achieved 
through a niobium pin, 2 mm in diameter, which could 
be moved perpendicularly to the plane of the billiards 
from outside the helium cryostat by a drive. The cou- 
pling strength is determined by the depths X\ and x 2 by 
which the niobium pin penetrates into the 71- and 72- 
stadium, respectively. For the strongest coupling, a sec- 
ond niobium pin, penetrating all the way through both 
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resonators, was added. Stronger coupling could have 
been obtained by using even more coupling pins. This 
was, however, not realized since it was the particular em- 
phasis of the present experiment to study the transition 
from the uncoupled case into the weakly coupled regime. 
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FIG. 2. The adjustable superconducting coupling. 
For each of the couplings we have measured the com- 
plete spectrum in steps of 10 kHz. In doing so, half of the 
microwave power was fed into each of the resonators. The 
states of the 71-stadium were not always visible through 
an antenna connected to the 72-stadium and vice versa. 
Therefore the spectrum of the entire system was con- 
structed by adding the spectra obtained through the 7 
antennas, four on the 71- and three on the 72-stadium, 
see Fig. |l[ 

A small frequency range of spectra at various couplings 
is shown in Fig. |^. One recognizes that the resonances are 
shifted by statistically varying amounts. The observed 
frequency shifts result from the coupling of the two cav- 
ities and from the perturbation of the electromagnetic 
field by the pin. The latter effect was investigated sepa- 
rately by inserting (or not) the pin into the 72-stadium 
only. Under that condition when only the field is per- 
turbed but no coupling is achieved all measured spectra 
(sec e.g. top of Fig. 4) displayed a clean 2-GOE behavior. 
The mean level spacing is D = 10.7 MHz. The resonance 
widths are of the order of 1 to 100 kHz. The high Q of 
10 5 — 10 6 together with the very good signal to noise ratio 
of up to 50 dB of the superconducting setup was obvi- 
ously necessary to detect the partly small shifts and thus 
the dependence of the level statistics on the coupling. 

We now turn to the analysis of the spectra. From 
the ansatz (^), the nearest neighbor spacing distribu- 
tion (NND), the S 2 -statistic and the A3-statistic can be 
obtained numerically, see and results below. Al- 

though the coupling parameter a/D can be determined 
by comparing numerical simulations with the data, it is 
convenient to have analytical expressions of £ 2 etc. as 
functions of a/D. French et al. |l^] and Leitner et al. 
Jl5[ have derived them for small coupling parameters us- 



ing perturbation theory. 

The analysis of the present data has been based on the 
£ 2 -statistic or number variance 



z{L) = ((n{L)-Lf). 



(2) 



Here, n(L) is the number of eigenvalues in an interval of 
length L. To obtain z(L), we divided the entire unfolded 
spectrum of length ND into Nl — N/L adjacent inter- 
vals of length L and took the average () over these. By 
looking at the correlation between z(L) and z(L + e), we 
convinced ourselves that z(L) and z(L + e) were statis- 
tically independent for e > 0.025 - at least in the range 
1 < L < 5. Calculating z(L) in that range in steps of 
e = 0.025 has provided M = 161 experimental numbers 
z(Lk), k = 1...M, that were statistically independent 
as is needed for the fit-procedure described below. The 
upper limit L < 5 is defined by the saturation of the 
E 2 -statistic |]l6| : Random matrix theory is known to 
describe spectral fluctuations of chaotic systems up to 
a maximum L max which is related to the length of the 
shortest periodic orbit. This sets L < 5 here. To check 
the influence of the saturation, we restricted the extrac- 
tion of a/D to L < 3. The slight change of the results 
was well within the errors. The lower limit of 1 < L is 
explained in the sequel. 
-10 
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FIG. 3. A small frequency range of three spectra with dif- 
ferent coupling. The arrows are intended to help to recognize 
the shifts of a few resonanc es. 

The expectation value z(L) is called £ 2 (L,A), where 
A = (a/D) 2 . According to Ref. jL5|, this function is 



£ 2 (L,A) = £ 2 (L,oo) 
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Here, £ 2 (L,oo) is the limiting function for the 1-GOE 
system. The parameter r is related to the ratio of dimen- 
sions in the GOE blocks of eq. (0). One finds r = 0.74 
in the present case. 

In order to estimate a/D, one has to know the prob- 
ability distribution Wk(z(Lk)\A), k = 1...M, of every 
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data point. By applying the "bootstrap method" to 
the set of Ni, k intervals from which z(L%) was calculated, 
we found w k to be a x 2 -distribution with average value 
as given by (||) and with roughly Ni, k degrees of freedom 
- which is reasonable. However, this was true only for 
Lk > 1. For L k < 1, no analytical representation of the 
distribution of z(L k ) was found. At the same time, the 
information on A is lost: for small L, the relative change 
of S 2 with A is of the order of L while the relative rms 
deviation of z is of the order of L 1 / 2 . Therefore, the 
analysis was restricted to L k > 1. 
1.5 




FIG. 4. The E 2 -statistic for three different couplings. The 
dotted line gives the 1-GOE and the dashed line the 2-GOE 
behavior. The solid lines result from the estimation of a/D 
described in the text and correspond to the first, the third 
and the last entry of Table N . 

The joint distribution TT(z|A) = J]fe w k (z(L k )\A) of 
all the z(L k ) was converted into the distribution W(A|;z) 
of A with the help of Bayes' theorem 



W(A\z) 



W(z\AMA) 



J dAW(z\A)n(A)' 
The a priori distribution of A was defined as 
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since this expression ensures that - at least for sufficiently 
large M - the entropy 
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W(z\A) \ 
f dAW(z\A)n(A) J 
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of W(z|A) is independent of A, whence A cannot be es- 
timated by a maximum entropy argument without any 
experiment. 

The center and the rms deviation of the distribution 
(El) were determined if it was Gaussian. This defines the 



best estimate and the error of A as well as a given in 
Table B. For three cases (part of) the data and the fit 
function (||) that results from the best estimate are given 
in Fig. (4). In the case of zero coupling, P(A\z) was not 
a Gaussian and the result given in Table | is an upper 
limit for the confidence of 68%. 

A % 2 -test - generalized to the case of non-Gaussian 
distributions w k - has shown that for the fits reported 
in Table [jj the analytical model (^) of [El| is compati- 
ble with the data. To check that the perturbation result 
(^) does apply to our case, we also performed numeri- 
cal simulations of the full model (|l|) using the procedure 
described in [jllj. Thereby we obtained values for a/D 
which cannot be affected by limitations of the perturba- 
tive calculation of jL5|. However, reassuringly, the results 
of both analyses are consistent within the errors. 

From the coupling parameter and the mean level spac- 
ing D — 10.7 MHz, one obtains the root mean square 
coupling matrix element a of the model ([l]). It corre- 
sponds to the rms Coulomb matrix element yj {H%) that 
is responsible for isospin mixing [ pT| . Figure |] shows 
that starting from 2-GOE behavior in the uncoupled 
case (a/D < 0.024) one moves through the weakly cou- 
pled case (a/D = 0.13) towards 1-GOE behavior. The 
strongest coupling (a/D = 0.20) realized here causes, 
however, still a relatively weak symmetry breaking of 
about the same size as the isospin symmetry breaking 
in 26 Al. The spreading width T L /D = 2ix(a/D) 2 which 
is a measure of how much the states of the two symmetry 
classes are mixed into each other, is also given in Table §. 
In the case of the strongest coupling e.g. one sees that 
a state of class 1 carries about 25% admixture of class 2 
and vice versa. This is the reason for the shifts observed 
in Fig. | 

The fact that the level statistics depend on a proves 
that the coupling block V in the Hamiltonian ([]]) is essen- 
tially filled with statistically independent elements - as 
we have assumed. If - on the contrary - V were a separa- 
ble interaction coupling only one specific configuration of 
the first symmetry class to one in the second class, then 
the level statistics would not change as a function of a. 

TABLE I. Mixing parameters for six different coupling 
strengths resulting from the Bayesian analysis described in 
the text. The penetration depths of the coupling pin into 
the resonators is given by (xi, X2) in mm. The results of the 
bottom row were obtained using a second coupling pin. 



Physical coupling 


a/D 


a (MHz) 




(0,8) 


< 0.029 


< 0.31 


< 0.0054 


(5,3) 


0.105 ±0.008 


1.12 ±0.08 


0.07 ±0.01 


(4,4) 


0.130 ±0.007 


1.39 ±0.07 


0.11 ± 0.01 


(5,8) 


0.173 ±0.006 


1.85 ± 0.06 


0.19 ±0.01 


(6,8) 


0.180 ±0.006 


1.93 ±0.06 


0.20 ±0.01 


(6,8) 


0.200 ± 0.006 


2.14 ±0.06 


0.25 ±0.01 
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FIG. 5. Fourier transforms of the spectra of the uncou- 
pled and a coupled system. The lengths of the periodic or- 
bits a, f3, 7 (shown in the insets) are indicated by the arrows . 
The lengths of the whispering gallery orbits are located within 
the hatched area. Orbit 7 touches the coupling pin which is 
marked by circles in the insets. 

We finally note that there is an experiment similar to 
the present one performed with elastomechanical reso- 
nances in quartz blocks [fLsf . Both experiments in prin- 
ciple allow semiclassical interpretations. The periodic 
orbits for fi~8|| are, however, very complicated. For the 
present relatively simple system, we have tried to detect 
periodic orbits that run back and forth through both of 
the coupled billiards. Systematic variations of the in- 
tegrated level density - including those caused by the 
bouncing ball orbits - have been removed as described in 
f|. The Fourier transform gf luc (x) of the remaining fluc- 
tuating part is expected to display the lengths x of the 
periodic classical orbits of the system 19 . In gf luc (x) 
obtained from the coupled stadia we have not been able 
to identify a peak introduced by the coupling and corre- 
sponding to an orbit running through both stadia. A very 
small part of the results is given in Fig. ^. In the range 
of x which is displayed, there are periodic orbits only in 
the 71-stadium. The shortest orbit of the 72-stadium is 
at 1.19 m. Introduction of the coupling changes gf luc (x) 
at every x - whether or not the orbits of length ss x 
came close to the coupling pin. This is expected be- 
cause \gf luc (x)\ 2 obeys a sum rule : The total intensity 
J dx\g' luc (x)\ 2 is given by the number of states. It is, 



however, interesting to see that a rather drastic change 
occurs in the vicinity of orbit 7. The peak at x « 1.03m 
(uncoupled case) splits into two (coupled case) and the 
interference minimum occurs at 7, i.e. the orbit 7 touches 
the coupling pin. 

In summary, the dependence of the spectral statis- 
tics on the coupling between levels belonging to different 
symmetry classes has been demonstrated for a system 
that simulates quantum chaos. Even subtle changes of 
the level statistics induced by small coupling parameters 
could be observed. The present experiment models mix- 
ing between any two symmetry classes. 
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